Gravitational clustering in a .D-dimensional Universe. 

T. Padmanabhan 1 * & Nissim Kanckar 2 ^ 

1 Inter-University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune 411 007, India 

2 National Centre for Radio Astrophysics, Post Bag 3, Ganeshkhind, Pune 411 007, India 



OS 



o 
O 



> 

in 

m . 
O- 
o- 

ov 

ON' 

Or 

o- 

a: 



X 



We consider the problem of gravitational clustering in a 
D-dimensional expanding Universe and derive scaling rela- 
tions connecting the exact mean two-point correlation func- 
tion with the linear mean correlation function, in the quasi- 
linear and non-linear regimes, using the standard paradigms 
of scale-invariant radial collapse and stable clustering. We 
show that the existence of scaling laws is a generic feature 
of gravitational clustering in an expanding background, in all 
dimensions except D = 2 and comment on the special nature 
of the 2-dimensional case. The D-dimensional scaling laws 
derived here reduce, in the 3-dimensional case, to scaling re- 
lations obtained earlier from iV-body simulations. Finally, we 
consider the case of clustering of 2-dimensional particles in 
a 2-73 expanding background, governed by a force —GM/R, 
and show that the correlation function does not grow (to first 
order) until much after the recollapse of any shell. 

I. Introduction 

The temporal evolution of a system consisting of a large 
number of particles interacting with each other via New- 
tonian gravity is a formidable problem to tackle. Such 
a system has, in fact, no "final" state of thermodynamic 
equilibrium as there is no limit to its phase space vol- 
ume. The phase volume available can be continuously 
increased by the separation of particles into a collapsed 
core and a dispersed halo, with the core becoming more 
and more tightly bound and the halo particles moving to 
larger and larger distances. In fact, the only stable con- 
figuration of such a system consists of a tightly coupled 
binary, with all the other particles at infinite distance. 
The above situation changes drastically if the background 
of the system (i.e. the space in which the particles are 
moving) is itself expanding. In this case, the expansion 
provides a "stabilising" influence on the evolution and 
can result in the formation of stable structures in which 
the effects of gravitational collapse are balanced by the 
expansion. Further, there also exists the possibility that 
the outmoving halos may be captured by other compact 
cores, leading to the build up of larger structures. Thus, 
the issue of gravitational clustering of a large number of 
particles appears more tractable in the case of a back- 
ground expanding (in general) with a time-dependant 
scale factor. This problem is, of course, of immense phys- 
ical relevance as there currently exists strong evidence 
that the matter density in the Universe is dominated by 
collisionless dark matter particles, which interact solely 
by gravity. Thus, if the length scales of interest are much 
smaller than the Hubble scale (and particle velocities are 
non-relativistic) , the formation of large-scale structure in 



the Universe is well described by the above picture, mak- 
ing it worthy of investigation. 

In the linear regime, where deviations from uniformity 
are small, perturbative techniques are used to obtain the 
time evolution of the system parameters, for example, 
the correlation functions. These methods, of course, fail 
in the quasi-linear and non-linear regimes and there is, 
as yet, no clear analytical picture of the behaviour of the 
system in these stages. However, in the case of 3+1 di- 
mensions, there exists a set of scaling relations relating 
the linear and non- linear correlation functions ( 0, pi). 
These scaling laws are not particularly well understood 
but appear to be validated by numerical simulations ( 

§, §)• 

A better understanding of a physical problem is some- 
times attained by treating it in a more general manner 
in an arbitrary number of dimensions as one may then 
be able to separate the generic features of the problem 
(as in, issues arising from the nature of the interaction) 
from results which stem from its dimensionality. This 
has been seen earlier, for example, in the case of the 
Ising model. In the current letter, we address the issue 
of gravitational clustering in a D-dimensional, expand- 
ing Universe and attempt to derive the scaling relations, 
using the well-known paradigms of scale-invariant radial 
infall in the quasi-linear phase and stable clustering in 
the non-linear regime. 

II. Clustering in D-dimensions 

We consider the case of a D-dimensional Universe, ex- 
panding with a time-dependant scale factor a(t). Fur- 
ther, we use the maximally symmetric Robertson- Walker 
metric (in D + 1 dimensions) and specialise to flat space 
(k = 0). The equations governing the evolution of a(t) 
are then (see S) 
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where pb and p are, respectively, the density and pressure 
of the background Universe and k is the constant in the 
Einstein equations, which, in general, can depend on the 
dimension D. We note that equations (|l|) and (||) are ob- 
tained by imposing the constraints of homogeneity and 
isotropy on the Einstein equations in (D + l) dimensions; 



the situation is thus manifestly isotropic. Further, it can 
be seen that the behaviour in the D = 2 case is likely to 
be special, due to the presence of the (D — 2) coefficient in 
the second term of equation (g) . We emphasise that the 
present work considers the clustering of D-dimensional 
particles in a Z3-dimcnsional expanding Universe. Ear- 
lier studies of 2-D clustering in the literature (see, for 
example, |J7j) treated the clustering of infinite needles in a 
3-dimensional background; this will be commented upon 
later. 

Given an equation of state, p(pb), one can solve equations 
(|l|) and (g) for the evolution of a(t). For pressureless 
dust, with p = 0, this implies that pf, is given by 
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Next, we incorporate the effects of perturbations on the 
smooth background by considering the evolution of the 
above system starting from Gaussian initial conditions 
with an initial power spectrum Pj„(fc). The two-point 
correlation function, £(a, x), is defined as the Fourier 
transform of the power spectrum. We will, for conve- 
nience, work with the mean two-point correlation func- 
tion, £(a,x), defined by 



where I is some other length scale. (Note that x (1 + £) 
is proportional to the number of neighbours of a given 
particle, within a Z?-dimensional sphere of radius ax; the 
above equation expresses the conservation of pairs in this 
sphere M . ) When the evolution is linear at all relevant 
scales, £(a, x) <C 1 and x ~ I. However, in the non-linear 
regime, £(a, x) S> 1 at some scale x] clearly x <^i I. Thus, 
the behaviour of clustering at a scale x is determined by 
the transfer of power from a larger scale I along the char- 
acteristics defined by equation (Q). This suggests that 
one should try to relate the true £(a,x) to the linear cor- 
relation function £^(a, I), evaluated at a different scale. 
The paradigm of scale-invariant radial collapse (see M, 
B) will be used to carry out the above procedure, in 
the quasi-linear regime. Consider the evolution of a 
D-dimensional, spherically symmetric, overdense region, 
containing a mass M . In general, such a region will ini- 
tially expand with the background Universe until the ex- 
cess gravitational force due to its enclosed mass causes 
it to collapse back upon itself. Thus, the radius of the 
region will initially rise, reach a maximum and then de- 
crease. The equation of motion for the radius, R(t), of 
such a £>-dimensional, spherical region is ( [pi) 
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and attempt to relate the exact £(a, x) in the quasi-linear 
and non-linear regimes to the mean correlation function 
calculated from linear theory, at the same epoch. 
The equation for the conservation of pairs can be written 
in /^-dimensions as 
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where v(a, x) is the mean relative pair velocity at scale x 
and time t. In terms of £(a, x), this gives 
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Here, h(a,x) — —v/dx, X = In x and A = In a. The 
above equation can be further simplified by defining 
F = In [x D (l + £)], yielding 
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where we have replaced for the mass, M, in terms of 
Si, I and ti, i.e. the initial density contrast, shell radius 
and time, respectively. We note that the acceleration, 
(d 2 R/dt 2 ), is proportional to l/i? D_1 , in D dimensions. 
The energy integral for equation ffl) is 
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The above expression is, of course, not valid for D = 2; 
this case will be treated separately later. We evalu- 
ate the energy constant, E, by requiring that the ve- 
locity satisfies the unperturbed expansion at t = ti, i.e. 
dR/dT = IHi = 2l/Dti, where H t = 2/Dt t is the Hubble 
parameter. This gives 
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At turnaround, R = Rm and dR/dt — 0. Thus 



The characteristic curves of this equation, on which F is a 
constant, satisfy the condition In [a; £) (l + £)] = constant, 
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or, since Si <C 1, 
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In the quasi-linear regime, we expect clustering to take 
place in a region surrounding density peaks of the linear 
regime, i.e. around regions such as the spherical region 
considered above. Making the usual assumption that the 
typical density profile around such a peak is equal to the 
average profile around a mass point, we can write this 
density profile as 
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Thus, the initial density contrast, Si(l) oc £,l(1), in the 
initial epoch, when linear theory is valid. Since Rm = 

IS { , clearly Rm oc £ l . In the scale- invariant, 

radial collapse picture, each shell can be approximated 
as contributing an effective radius proportional to Rm- 
Taking the final effective radius, x, as proportional to 
Rm, the final mean correlation function is given by 



mean correlation function have the relation £ oc S 2 , in 
the linear regime; this arises from the definition of £ as 
the Fourier transform of S 2 . Now, in a flat D-dimensional 
matter-dominated Universe, S oc a D ~ 2 ( ||), in the lin- 
ear regime. Thus, £ L oc Sf oc a 2 ( D ~ 2 \ We can hence 
write £i(a, I) — a 2 ( D ~ 2 'Q \l D \ , where Q is a function of I 
alone. But, l D « x D l Nh (a,x) = x D a D F(ax) = r D F(r). 
This implies 
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Considering the above relation as a function of a at con- 
stant r, we obtain 
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where B and C are constants. One must hence have 
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Thus, in the extreme non-linear regime, we obtain 
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Thus, the final correlation function £q L at x is the 
D/(D — 2) power of the initial correlation function at 

I, where l D oc x D [£^(0] oc x d £,ql{x), which is 

the form required by equation (jg), if £q L 3> 1. 
Let us turn our attention, next, to the non-linear regime; 
here, we use the ansatz of stable clustering (h — ► 1 for 
S, ^> 1), which is physically well-motivated as it seems 
reasonable to expect stable, bound systems to form un- 
der the joint influence of gravity and the expansion ( 
P|). Such systems would neither expand nor contract 
and would hence have peculiar velocities equal and op- 
posite to the Hubble expansion, i.e. v 1 = —ax 1 . 
(We emphasise that the stable clustering hypothesis is 
an ansatz and might well fail if mergers of structures 
are important.) Stable clustering requires that virial- 
ized systems should maintain their densities and sizes 
in proper co-ordinates, r — ax. This would require the 
correlation function, in D dimensions, to have the form 
£ NL (a, x) — a D F(ax), where F is some function and the 
aP factor arises from the decrease in the background den- 
sity. We assume that £,^^(a,x) is a function of t; L (a,l), 
where l D « x D £, NL (a,x), from equation (||). This rela- 
tion can be written as 
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The above analysis shows that the exact mean correla- 
tion function can be expressed in terms of the linear mean 
correlation function by the relation 
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(quasi-linear) (23) 
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In the case of 3 dimensions, the above equations reduce to 

£(a, x) oc [Cl(°;0] an d [Ci( a >0] m ^ ne quasi-linear 
and non-linear regimes, respectively' these are in reason- 
able agreement with simulations ( 0], jqj; see, however, 
0). Thus, scaling laws appear to be a generic feature of 
gravitational clustering in an expanding background, in 
all dimensions (except D — 2). 

The 2-dimensional case is special and different from all 
other dimensions, as, in this case, equation (M) gives 



where U(p) is an unknown function of its argument, 
which is to be determined. The density contrast and 



d 2 R 

~aW 



= 



(25) 



i.e. the correlation function does not evolve at all in 2-D 
and no structures are formed. The special nature of the 
D = 2 case appears to arise from the structure of the 
Poisson equation, due to the fact that it contains deriva- 
tives of the second order. It has been noted earlier that 
scaling relations can arise in 2-D numerical simulations 
of gravitational clustering ( @). These simulations, how- 
ever, do not consider the case of 2-D particles evolving in 
a 2-D background but instead treat the system as con- 
sisting of a set of infinitely long, parallel needles, in a 
3-D expanding background, with particles considered as 
arising at the intersections of the needles with any plane 
orthogonal to them. The mass elements in the needles 
interact via the usual 1/r 2 force; however, the interac- 
tion between the needles themselves is governed by a 1/r 
force. The simulations consider the clustering of the nee- 
dles by taking a slice through a plane orthogonal to them. 
This situation is, of course, clearly anisotropic as cluster- 
ing is considered as taking place in 2 dimensions while 
the background Universe expands in 3 dimensions. In 
this case, one can show numerically ( |l]]) that Rm = I /$i 
in the quasi-linear regime and arrive at the relevant scal- 
ing relations by an analysis similar to the above. 
We finally consider the case of 2-D particles in a 2-D ex- 
panding background and assume a force law 
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The above form, of course, cannot be obtained in a self- 
consistent manner, by taking limits of the Einstein equa- 
tions for an expanding Universe; however, it has the form 
d 2 R/dt 2 oc i? _1 , the correct dependence for D = 2 in 
equation (m . The energy integral of equation (|2q) is 



Thus, the correlation function does not grow (to first or- 
der) until much after turnaround. This is an interesting 
result which needs to be verified by simulations. 
The above situation is similar to the case of a 3-D net- 
work of parallel cosmic strings (space around them is flat 
and pb ex a~ 2 ). Structure will clearly not grow if the net- 
work is expanding uniformly; however, if the strings had 
some initial peculiar velocities, then two strings which 
pass on either side of a third would move towards each 
other because the angle around the third string is less 
than 27r. 

III. Conclusions 

In the present paper, the standard paradigms of scale- 
invariant radial collapse and stable clustering have been 
used to derive scaling relations connecting the exact 
mean two-point correlation function and the linear mean 
correlation function, for a D-dimensional expanding Uni- 
verse. The existence of scaling laws is found to be a 
generic feature of gravitational clustering in an expand- 
ing background, in all dimensions, except D = 2. Fur- 
ther, the scaling laws derived here are in agreement, in 
the 3-dimensional case, with scaling relations obtained 
earlier from ./V-body simulations (£ ex ££ and £ ex £ L ' , 
in the quasi-linear and non- linear regimes, respectively). 
Finally, we have considered the case of clustering of 2-D 
particles in a 2-dimensional expanding background, gov- 
erned by a force -GM/R, and show, by a similar analysis, 
that the correlation function does not grow, to first or- 
der, until much after turnaround of any shell. 
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The mass M, enclosed within the shell, is given by 



M = (1 + 5i)p b nl 2 



(28) 



ting 



where I is the initial shell radius. Here, we will assume 
Pb ex t~ 2 , the usual result for D dimensions. Sett 
dR/dt — Hil initially, and using Hi ex !/£,, we obtain 



R M =lexp[-A(l + 5i)] 



where A is a constant; i.e. Rm/1 = const., to leading 
order in £,-. Again taking the final effective radius, of 
the 2-dimensional shell, as proportional to Rm, the final 
mean correlation function is given by 



(27) 
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